Report Number — NPS55-83-026 



FEDDOCS 

D 208.1 4/2 :NPS-55-83-026 



NAVAL 



POSTGRADUATE SCHOOL 

Monterey, California 




WHOLESALE PROVISIONING MODELS: 
MODEL DEVELOPMENT 



F. Russell Richards 
Alan w. McMasters 



September 1953 



Final Report tor Period January 1982- June 1983 

Approved for public release; distribution unlimited 

Prepared for: 

Commanding Of f icer 

Fleet Material Support Office 

Mechanics burg, PA 17053 



NAVAL POSTGRADUATE SCHOOL 
Monterey, California 



Rear Admiral J. J. Ekelund 
Superintendent 



David A. Schrady 
Provost 



The work reported herein was supported with funds provided by the Navy 
Fleet Material Support Office, Mechanicsburg , Pennsylvania. 

Reproduction of all or part of this report is authorized. 



This report was prepared by: 



UNCLASSIFIED 

SECURITY CLASSIFICATION OF THIS P ACE (Whan Dmtm Entarad) 



REPORT DOCUMENTATION PAGE 


READ INSTRUCTIONS 
BEFORE COMPLETING FORM 


1. REPORT NUMBER 

NPG55-83-026 


2. GOVT ACCESSION NO. 


3. RECIPIENT'S CATALOG NUMBER 


«. TITLE (end Subtitle) 

WHOLESALE PROVISIONING MODELS: 
MODEL DEVELOPMENT 


3. TYPE OF REPORT 6 PERIOD COVERED 

FINAL REPORT 

1 January 82-1 July 83 


S. PERFORMING ORG. REPORT NUMBER 


7. AUTHOR(«j 

F. Russell Richards 
Man W. McMasters 


• . CONTRACT OR GRANT NUMBSRf*; 


9. PERFORMING ORGANIZATION name and aodress 

Naval Postgraduate School 
Monterey, CA 93943 


10. PROGRAM ELEMENT, PROJECT, TASK 
AREA * WORK UNIT NUMBERS 

N0036783WRM4674 


11. CONTROLLING OFFICE NAME AND ADDRESS 

Commanding Officer, Fleet Material Support Office 
Mechanicsburg, PA 17055 


12. REPORT o’aTE 

September 1983 


IS. NUMBER or PACES 

60 


14. MONITORING AGENCY NAME A ADORESSf// dlffaeant Irom Controlling Oftl c*) 


is. SECURITY CLASS, (o! thia raport) 

Unclassified 


I3«. OeCL A$5IF| CATION/ DOWN GRADING 
SCHEDULE 



16. DISTRIBUTION STATEMENT (of thia Report) 



Approved for public release; distribution unlimited. 



17. DISTRIBUTION STATEMENT (of tha abatemct ontarad In Stock 20, tt dlffarant from Raport) 



18. supplementary notes 



19. KEY WORDS (Confinu* on tavaraa aid* it nacaaamry and identity by block numbar) 

Provisioning Mean Supply Response Time 

Availability 

Inventory Models 

2G. ABSTRACT (Conttnuo on revora a aid a If nacaaamry mid Idantlfy by block mtmbrnr) 

Alternative wholesale provisioning models for secondary items are developed^ 
The models allocate optimally a given provisioning budget with such objective 
functions as units short, supply material availability, mean supply response 
time and availability. The models allow the inclusion of essentiality 
weights to discriminate component criticality. A procedure is also discussed 
for determining the minimum budget required to satisfy a specified level of 
performance . 

— . — — ■ 



DD 1 jan M 73 1473 EDITION OF 1 NOV 85 13 08S0LETE 

S' N 0102- LF- 0U* 6601 



UNCLASSIFIED 



SECURITY CL ASSlFlCATlON OF THIS PAGE 'Whan Dmt a Jntmrad) 



WHOLESALE PROVISIONING MODELS 



MODEL DEVELOPMENT 



by 



F . Russell Richards 
and 

Alan W. McMasters 



NAVAL POSTGRADUATE SCHOOL 
MONTEREY, CALIFORNIA 



VC 
A & 3 
R ys 

i ?%3 



TABLE OF CONTENTS 



ABSTRACT 

1. INTRODUCTION 

2. THE CURRENT SPCC PROVISIONING MODEL 

Development of the Budget Constraint 

The Variable Threshold Model for Budget 
Implementation 

3. MODEL BACKGROUND 

Introduction 

The Provisioning Interval 

Provisioning Objective Functions 

Form of the Models 

Range Rules 

4. THE DEMAND DISTRIBUTION 

5. UNITS SHORT MODELS - SMA 

Introduction 

Model I - Units Short 

Model II ~ Modified Units Short 

Model III - Supply Material Availability 
Summary 

6. TIME-WEIGHTED UNITS SHORT 

Model IV - TWUS 

Model V - Mean Supply Response Time 

7. AVAILABILITY 

8. COMMENTS ON THE MODELS 

Introduction 

Random Protection Interval 

Theoretical Properties 

A Companion Problem as a Budgeting Tool 

9 . SUMMARY 

APPENDIX A: INITIAL DEMAND FORECASTING 

APPENDIX B: THE C OSD IF FORMULA 

The Philosophy 

A Critique 

APPENDIX C: AVAILABILITY vs, MSRT 

LIST OF REFERENCES 



i 

X 

5 

5 

6 

10 

10 

11 

12 

13 

14 

15 

13 

13 

13 

22 

23 

25 

27 

27 

31 

34 

33 

38 

38 

40 

41 

43 

45 

49 

49 

51 

5b 

58 



ABSTRACT 



In this report we develop alternative wholesale provisioning models for 
secondary items with an aim toward use by both the Ships Parts Control Center 
(SPCC) and the Aviation Supply Office (ASO) . General background into the 
wholesale provisioning problem is given, and the present provisioning models 
in use at SPCC are summarized. We review the demand forecasting problem, and 
we address the issue of demand distribution assumptions. Then various 
alternative wholesale provisioning models are developed. These models 
attempt to allocate optimally a given provisioning budget with such objective 
functions as (1) units short, (2) essentiality-weighted units short, (3) 
essentiality-weighted requisitions short, (4) supply material availability 
(SMA), (5) time-weighted units short, (6) essentiality-weighted mean supply 
response time, and (7) system availability. We conclude with theoretical 
comments about the alternative models and discussion of a companion problem 
that can be solved to determine the minimum budget required to attain a 
specified level of performance. 

This is the first of three reports on the wholesale provisioning problem. 
The second report will present results of a numerical evaluation of the 
proposed alternative models and the current model. The thira report addresses 
the issue of phased provisioning and offers a heuristic method for determining 
how to phase budget expenditures over time. 



1 . INTRODUCTION 



In the spring of 1982, the Naval Supply Systems Command (NAVSUP) asked 
the Naval Postgraduate School to develop improvements to the existing 
peacetime wholesale provisioning models for secondary items used by the Ships 
Parts Control Center (SPCC) and the Aviation Supply Office (ASO) . This effort was 
motivated by NAVSUP T s Resolicitation Project, the major objective of which 
is to acquire new computer hardware. However, it also provides an opportunity 
to take a hard look at existing models in use by NAVSUP f s Inventory Control 
Points (ICP) at SPCC and ASO and to make changes as desired before the ICP 
software is installed in the new computers. These changes should include 
wholesale and retail model improvements, both in provisioning and 
replenishment of secondary items . The most important questions to be 
answered by the Naval Postgraduate School are: 

1. Do the existing models provide the best supply support for the Navy’s 
weapon systems? 

2. How should those models which do not provide the best support be 
changed? 

In responding to these questions, it makes sense to begin with "the 
beginning" of supply support — provisioning for a weapon system — because, once 
an improved provisioning model is developed, the replenishment models can be 
modified in such a way as to provide a smooth transition to continued 
wholesale supply support. This report will therefore concentrate on the 
provisioning models. 

Wholesale supply support provides backup or "systems" stock of repairable 
and consumable items for the retail levels. At present the retail levels use 
models such as COSALs (Coordinated Shipboard Allowance List) and A7CALS 
(Aviation Coordinated Allowance List), These are stocks of spare secondary 
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items (items which may fail in use) designed to provide sufficient support for 
a ship deployment period of, say, 90 days. Whenever an item is used by a ship 
or a squadron, a requisition is immediately submitted to the wholesale system 
for a replacement. Thus, the COSALs and AVCALs serve only as emergency 
protection. All demands ultimately must be satisfied by the wholesale stocks. 

The wholesale provisioning model provides protection until the wholesale 
replenishment model can initiate the first replenishment buy (or repair 
quantity for those secondary items which can be repaired) and that buy is 
manufactured and delivered. Thus, wholesale provisioning buy protection 
begins on the date of preliminary operational capability (POC) and continues 
for at least a replenishment procurement (or repair) leadtime. 

Actually, it is expected that a wholesale replenishment buy will be 
initiated some substantial time after POC. This will occur when the reorder 
point established by the replenishment model is reached. The initial reorder 
point is determined differently by the two ICP's and it may be based on the 
initial expected demand which is predicted as part of the provisioning 
process. However, after six months, the reorder point begins to reflect the 
actual demands placed on the supply system. This may result in the value of 
the reorder point decreasing for a year or two. 

The provisioning buy must be made early enough so that the stock will be 
on hand by the POC date or shortly thereafter. Thus it must be made at least 
a procurement leadtime prior to that date. Typical values of this initial 
leadtime appear to be between 2.0 and 2.5 years. With such long leadtiraes the 
quality of forecasted installation schedules of a new system may be poor. 

And, when protection must be provided for a subsequent time which is at least 
as long as a replenishment leadtime, the quality of the forecasted demands 
over that time can also be poor. 
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Provisioning guidance is provided by the Department of Defense 
Instruction (DODINST) 4140.42 [7]. The authors of [7] were of the opinion 
that much stock of secondary items was procured which was not subsequently 
needed. Thus, in an attempt to remedy this perceived situation, a model was 
developed which would provide only a very conservative quantity of wholesale 
stock for those items which were new to the supply system. For items already 
stocked in the system, they allowed no additional stock to be procured 
specifically for a new system. However, if the anticipated increase in demand 
was sufficiently large, it could be entered into the demand forecasting 
process by the inventory manager and its influence could trigger an earlier 
replenishment buy than would have occurred without it. 

The authors of [7] did appreciate that improved models for provisioning 
new items would evolve and stated that: 

Inventory optimization models may be used as a basis for stockage in lieu 
of the model in this Instruction if the following criteria are met: 

(1) An optimization technique, developed to minimize system downtime or 
time-weighted requisitions short is used. (2) No lower limits are placed 
on requirements that will result in stockage as a cemand-based item, of 
an item that would not be stocked without the lower limits. (3) A 
monetary base point is determined indicating the approximate value of the 
requirement in accordance with the instructions contained in enclosure 2. 
The second criterion may be omitted if only those items meeting the 
stockage criteria in accordance with this Instruction are considered in 
the optimization model. The longest applicable Program F orecas t Period 
may be used in the optimization model, after the monetary base has been 
es tablishea. 
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Such improvements have indeed been developed by NAVSUP and are in use at 
both SPCC and ASO . However, as was stated above, NAVSUP is interested in 
developing even better models. This report represents the first step and will 
concentrate on a model to improve the initial supply support for new secondary 
items that will be needed. The issue of how to handle items already stocked 
will be the subject of a later report. 

This report begins with a discussion of the existing wholesale 
provisioning model used by SPCC since it is more mathematically sophisticated 
than the ASO model and it is also easier to understand. Next, the demand 
process is examined and several alternative models intended for use by both 
SPCC and ASO are developed. These alternative models develop a budget using 
the COSDIF process and then allocate that budget so as to optimize system 
perf ormance with respect to such measures as supply material availability, 
mean supply response time, and availability. A discussion of some of the 
theoretical properties of the models is provided and recommendations based on 
these properties are made. 

Three appendices are included. The first presents the details of the 
forecasting models from [7] and from the SPCC and ASO models. The second 
discusses the cost difference or COSDIF formula which is required by [7] for 
developing the budget constraint. The last examines a widely held 
misconception about the relationship between mean supply response time and 
availability. 

This is the first of three reports cn wholesale provisioning. The second 
report describes the computational algorithms used for provisioning with each 
of the alternative models and provides some numerical evaluations of the 
models. The third report examines the issue of phased provisioning and offers 
a heuristic procedure for determining how to divide a provisioning budget over 
1 period of time. 



4 



2. THE CURRENT SPCC PROVISIONING MODEL 



DEVELOPMENT OF THE BUDGET CONSTRAINT 

The model that SPCC uses to determine initial wholesale stocks for new 
items is a variant of the model proposed in reference [7]. The application/ 
operation which implements the model is designated as D55 [19]. The model 
makes range and depth decisions separately for each item in the system for the 
purpose of developing a budget constraint. 

The initial phase of the process is to develop a budget constraint. The 
first step of this phase is the determination of a range of items to be 
considered for stockage. This is done by computing the COSDIF value for each 
item. The purpose of the computation is to decide if an item is to be stocked 
as a demand-based item. The COSDIF value is obtained from a formula which was 
stated in reference [7]. The current version of COSDIF used by SPCC is 
presented in reference [19] and is reproduced in Appendix B of this report. 
Appendix B also presents the motivation for the formula and a critique of the 
procedure. Elements of the COSDIF include the expected costs of stocking the 
item and having no demands over two years, the expected costs of stocking an 
item and experiencing demand for one year, and the expected costs of not 
stocking an item for one year and having to make spot buys when demands occur. 
The difference between the sum of the first two elements and the last element 
is the "cost difference" or COSDIF. If the value of the COSDIF is negative 
then the item will be stocked; otherwise it will not be stocked insofar as 
budget development is concerned. 

If an item has been determined by the COSDIF as appropriate to stock, 
its depth of stockage is next computed. The depth is taken to be equal to the 



demand expected over a period of time consisting of the forecasted replenish- 
ment procurement ieadtime plus one quarter. The extra quarter of demand is 
viewed as safety stock. The procedure used to develop the initial demand 
prediction is discussed in Appendix A. 

Those items which fail to satisfy the COSDIF criterion are next 
re-examined to determine if they are coded as insurance or numeric 
stockage objective (NSO) items. If so, the depth of each item is taken to be 
one minimum replaceable unit (MRU). 

The cost of buying the amount of each item to the depth specified above 
is then determined, and the total dollar value of a provisioning package is 
obtained by summing these procurement costs over all of the NSO items, the 
insurance items, and those items passing the COSDIF test. This sum is then 
established as "the budget constraint" in accordance with reference [7]. The 
dollar value so determined is considered as a firm upper limit on the amount 
of money that can be spent to purchase those items that are to be stocked. 

The actual range and depth of items that must be stocked, however, may 
deviate from the values used to determine the budget. As was noted in chapter 
I, an alternative provisioning model may be used to implement the range and 
depth provided the budget constraint as determined above is satisfied and 
the model objective is consistent with minimization of system down time or 
time-weighted requisitions short. The determination of the budget constraint 
is therefore an important element. 

THE VARIABLE THRESHOLD MODEL FOR BUDGET IMPLEMENTATION 

SPCC has developed a model, called the Variable Threshold Model, for 
making the budget implementation. In that model, all new items are considered 



to be candidates for stockage and are re-evaluated based on unit cost and the 
expected demand during the procurement leadtime. 

The first step of the model is to compute the variable threshold value of 
each item. This value is computed using the variable threshold formula: 

V(i) = (l - exp(-Di)] / C. 

where is the unit cost of the item and D-j_ is the expected demand during the 
procurement leadtime (see Appendix A). The variable threshold value is the 
mechanism used by SPCC to rank the items in a provisioning package according 
to desirability of stockage. Those items having a high value are placed at 
the top of the list. The numerator of the threshold value is the probability 
of one or more demands during a procurement leadtime when the demand is 
Poisson distributed. It comes from an analysis by Burton and Jaquette [2] in 
which the term represented the decrease in the expected number of backorders 
at a random time resulting from an increase in units stocked from zero to one 
for an item that has Poisson distributed demands. Reference ( 2 j addressed 
only very low-demand items which are currently viewed as insurance or NSC 
items . 

The variable threshold formula is a typical "benef it-to-cos t " ratio used 
in the first step of a marginal analysis allocation scheme. Similar 
expressions are found in the alternative models which we present later. 

The next step is to compute the depth for each item based on a risk 
formula similar to that used in the replenishment models. The formula is 

RISK = IC / (IC -r X S E) 

where IC is the annual holding cost, \ s is the shortage cost, and E is the 
item essentiality. 
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Let H(x) represent the probability that the leadtime demand for a given 
item is greater than or equal to x. Then, the depth for the item is taken to 
be the smallest value of x such that H(x) is not greater than the calculated 
RISK value. Thus, if x^* represents the depth for item i, 

x . * = min x; 

1 

1C. 

such that H . (x) < RISK = — — -- v ■ , 
i — IC. + A.E. 

i ii 

The actual probability distribution used for the calculation depends on 
the expected annual demand frequency. The Poisson distribution is used if the 
expected annual demand is less than or equal to one; the normal distribution 
is used when the expected annual demand or procurement leadtime demand is 20 
or more; and the negative binomial distribution is used when the expected 
annual demand is between one and 20. The actual depth is constrained to be 
no more than two years of expected demand if the item is consumable and no 
more than the expected f, net M attrition demand during a procurement leadtime 
plus one quarter if the item is repairable (see Appendix A) . 

The final step is to allocate the budget constraint among the items. The 
part of the total budget required to purchase the MRU's for the insurance and 
NSO items is reserved for these items and is subtracted from the total budget 
available. The remaining dollars are next allocated by successively 
purchasing the depths determined above for the items ranked in order of 
highest to lowest variable threshold value. The item with the highest 
variable threshold value is purchased in the depth specified above; the 
remaining budget is decremented; the item with the second highest variable 
threshold value is purchased next. This procedure continues until the budget 
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is depleted or until every item is purchased at the specified depth. Any 
items already designated as insurance or NSO items that qualify for stockage 
under the above procedure will first have the depth quantities reduced by the 
already funded MRU f s. 

A comparison of the SPCC budget allocation approach with that described 
in reference [7J was done by the Fleet Material Support Office (FMSO) and 
reported in reference [9]. The results of the study showed that the variable 
threshold model provided more range and less depth than did the model of 
reference [7] for the same budget amount. 
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3. MODEL BACKGROUND 



INTRODUCTION 

The provisioning model currently in use at SPCC was discussed in the 
previous chapter. That model is a variant of the model proposed in DODI 
4140.42. The provisioning procedure described in DODI 4140.42 is an aggregate 
of the Army's replenishment range model (COSDIF) and an attempt to 
conservatively forecast leadtime demand (under the assumptions that the 
components have reliabilities equal to their design values). That instruction 
provides uniform guidelines for the military services under which the services 
could develop their own provisioning models for carrying out the budget 
implementation The instruction does, however, allow only certain types of 
models. Paragraph IV. D of the instruction states: 

“Mathematical models may be used which may provide for a different mix of 
inventory than that stated herein, providing a financial base is 
established in accordance with the policies outlined in this instruction 
and an objective to minimize system downtime or time-weighted 
requisitions short is included in model.” 

we examine models in this report which focus on these types of objectives. In 
particular, we develop models which optimize with respect to: 

1 . expected number of stockouts 

2. expected time-weighted units (requisitions) short 

3. mean supply response time, and 

4 . availability. 

Before presenting these models, it is appropriate to examine some background 
issues «7hich affect all of the models that follow. 
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THE PROVISIONING INTERVAL 



There are four key dates for consideration in the provisioning problem: 

1. POC: the date of planned operational capability (part of the 

provisioning package should be in place by this date) , 

2. POC - PCLT: a procurement leadtime prior to POC , 

3. POC + TR: the time beyond POC at which the first replenishment buy 
is made , and 

4. POC + TR + PCLT: the time at which the first replenishment buy is 

received. 

The primary concern of provisioning is to determine the number of units 
of each component of a system to be purchased at time POC - PCLT (and possibly 
purchase at additional times between POC - PCLT and POC). In order to reduce 
the likelihood of stockouts * the provisioning buy made at time POC - PCLT 
should be sufficiently large to satisfy anticipated demands in the interval 
(POC, POC + TR + PCLT). If good forecasts of demand during that interval and 
of the times TR and PCLT were available, the solution of the provisioning 
problem would be fairly routine. However, the provisioning problem is 
characterized by a great deal of uncertainty with respect to the failure rates 
(3RF f S) of new equipments, the procurement leadtime, the time TR, and even the 
installation schedule of new equipments. 

During the time interval from POC to POC ~ TR + PCLT, the population of 
equipments will grow as additional units are installed in the operational 
community. This will cause the aggregate failure rates to increase over the 
time interval of concern in the provisioning problem. Actual installation 
schedules are subject to a great deal of uncertainty; the 3RF estimates often 
represent only guesses at the failure rates since there are little or no 
operational data from which to estimate the failure rates; and the time lags 
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are very large (the time interval from POC - PCLT to POC -r TR + PLT may be on 
the order of four years or more). Thus it is not surprising that the 
provisioned quantities frequently don f t match the demands very well. Faced 
with such uncertainties, no analytical model will cure all of the provisioning 
ills. However, the use of rational provisioning procedures should minimize 
such problems and, in addition, lead to improved replenishment models. 

PROVISIONING OBJECTIVE FUNCTIONS 

A key ingredient of DODI 4140.42 is the budget constraint for the 
procurement of items to be provisioned. No flexibility exists in the method 
used to determine the budget. (The actual values of the parameters used in 
the COSDIF formula can, however, be selected by the users.) As was noted at 
the start of this chapter DODI 4140.42 does allow limited flexibility in the 
type of provisioning model used to allocate the given budget among new items 
being considered. It can either minimize time-weighted units short or 
minimize system downtime. However, models such as the variable threshold which 
actually do neither have been approved. As a consequence, many measures of 
system performance have been and continue to be considered for objectives of 
provisioning. Some of the more popular are availability, probability of a 
stockout, fill rate, backorder accumulation, system downtime, time-weighted 
requisitions short, essentiality-weighted shortages, mean supply response 
time, and inventory costs. The use of costs as the principle measure of 
performance is complicated by the impossible task of quantifying such elements 
as stockout costs and obsolescence costs. The use of true system availability 
models for the wholesale level is complicated be the need for information at 
the operational level on a systera-by-system basis. That information is 
currently not available. Such models are therefore not presently feasible for 
the wholesale level. 
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Some implementations of availability models for stockage decisions define 
availability for individual components (as opposed to weapon systems) in terras 
of the ratio of mean time between failures (MTBF ) to the sum of MTBF , mean 
time to repair (MTTR) , and mean supply response time (MSRT) (or other similar 
expressions); see, for example, references [ 3] and [23] . The proponents of 
such availability models argue that MSRT is the only terra which pertains 
directly to the provisioning stockage decision. Consequently, these 
availability models really attempt to minimize MSRT. (Appendix C provides a 
simple example to show that minimization of MSRT does not necessarily yield 
the same solution as maximization of system availability.) 



FORM OF THE MODELS 

In all of the models which we develop, the optimization problem is of 
the form: 

minimize (maximize) T f.(s.) 

.11 

i 



subject to 



l c . s . \ 

. 11 — 



< B, 



where f^(si) is the performance level for item i when s^ units are stocked, c^ 
is the unit price of item i, and 3 is the total amount budgeted for the 
provisioning package. The form of the optimization problem implicitly assumes 
separability; that is, the total system performance can be separated into a 
function which is the sum of the performance generated by each of its 
components. This form may not always be appropriate because it ignores the 
system configuration and any interaction between components. However, we feel 
that such an assumption is necessary to assure mathematical tractibility and 
implementation feasibility. It is also in keeping with the piece-parts 
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support view which is traditional in the supply system. True "systems 
support" models are still a long time away. 

When the performance functions are separable and well behaved the 
mathematical solution to the optimization problem can readily be obtained 
using such straightforward techniques as the generalized Lagrange multiplier 
approach of Everett [ 3] or marginal analysis [12]. 

RANGE RULES 

Ail of the models which we will address solve simultaneously the range 
and depth problems. There is no explicit: range rule. Instead, zero depth is 
equivalent to a no-stockage decision. Insurance items and numeric stockage 
objective (NSO) items will have to be handled separately from demand based 
items just as they are today. For the items which are demand based the models 
require an assumption about the distribution of demand. The next chapter 
supports the use of the Poisson probability distribution with an expected 
demand rate which depends on the BRF estimate and a time-weighted installed 
populat ion ( TWAMP ) . 
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4 . THE DEMAND DISTRIBUTION 



One of che most critical elements in any stockage model is the estimation 
of the total number of failures (demands) for a component which will be 
generated during the provisioning protection interval. The proper way to 
handle demand uncertainty is to model the demand process as a random process 
with a probability law from a known family of probability distributions and 
having parameters which are estimated from the specific information that is 
available for each item. In the provisioning problem, this information is 
pretty skimpy, consisting essentially of contractor provided estimates of 
failure rates (TRF ? s) and the schedule of installations. 

This chapter discusses the probability distribution that should be used 
for the provisioning stockage decisions. We show that, under the assumption 
of a constant failure rate, the total demand for a given item over any time 
interval is Poisson distributed with parameter X which depends on the 
time-weighted average month's program (TWAMP). 

Consider a single component with constant failure rate a during the 
interval ( 0 ,T ] • Let N be the total number of installations of the component 
during ( 0 ,T ] , and let 0 = t \ < t2 < . . . < t^j be the installation times of the 
N units of the component. Let X be the random number of failures of unit i 
in (0,T]. Then, because of the constant failure rate, the distribution of 
is Poisson with parameter a(T-tj_) (see, for example, Ross [21]); i.e., 

P(Xi = x) = (a(T-ti) ) x exp(-a(T-tj_ ) ) / x! x=l,2,... 

Let us now consider the total number of demands generated by the i ^ and 
j tn installed units. If we assume that the failures of the i c ^ and units 

are independent, the probability distribution of the sum of Xj_ and Xj is given 
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by the convolution of the distributions of and X j . Thus, 



m 

P(X. + X. = m) = l P(X. = m-k)P( X. = k) 

1 J k=0 3 1 

m m-k k 

^(a(T-tj)) exp(-a(T-tj ) ) (a(T-ti)) exp(-a(T-ti) ) 

k= ° (m-k)! kT 

V 

If we let p = (T-t^) / (2T-ti-tj), then 1-p = (T-tj) / (2T-ti-tj) and the 
above summation is recognized to be the sum of binomial probability masses. 
Hence, that sum is unity, and the distribution of the sum of and Xj is 
Poisson with parameter a( 2T-t ^-t j ) . Let D = X^ + X 2 + • . . + Xjj be the total 
demand over all installed units during the interval (0,Tj. The above result 
is easily extended to show that D is Poisson distributed with parameter 
m 

a (NT- T t-^). Let f(t) be the total number of installations at time t. Then 
i=l 

f(t) is a step function with jumps at times t \ , t2 , . t^ • The integral of 

N N 

f(t) over the interval is therefore ^ (T-t^) = NT - 5 t^. But this integral 

i=l i=l 

when normalized by T is the definition of the TWA11P . (TWAMP as defined in 
\ppendix A treats time as discrete in units of months and therefore must make 
certain adjustments to accomodate installations occuring during a month. The 





16 



differences are minor). Thus, the total demand over the program time base 
(0 ,T] is Poisson distributed with parameter aT*TWAMP. Figure 4.1 illustrates 
this for the case in which 9 total installations occur over the interval 
(0,T] . 

If all of the installations took place at POC (t[ = t2 = . . . = tjj = 0) , 

the process describing the total demand is a homogeneous Poisson process with 

rate X = N a. However, because of the phased installations, the total demand 

process has a mean rate which is a nondecreasing function of time; i.e. , the 

demand process is a non-homogeneous Poisson process. The mean rate of demand 

t 

at time t is A(t) = a f f(t)dt. 

0 

In the next chapters we develop alternative models for determining the 
amount of stock to be provisioned for each component of a provisioning 
package. In all of those models we assume a Poisson distribution for the 
total demand for each item. The parameters \ ^ referred to in those models 
will be assumed to be the "TWAMP'ed" values discussed above. 




Figure 4.1: Time-Weighted Installations 



17 



5. UNITS SHORT MODELS (SMA) 



INTRODUCTION 

Consider a weapon system composed of n components. Let s^,s^,...,s^ 
be the number of units of components 1 through n, respectively, to be 
provisioned. Let c_^ be the unit cost of component i and let 3 be the 
total budget available to fund spares for the provisioning package. (B will 
be determined by the COSDIF process.) As discussed earlier, the provisioning 
problem is interested in the time interval (POC, POC + TR + PCLT) . However, 
for notational convenience we will refer to the interval as (0,T ]. Let 
X.(T.) be a random variable describing the total number of demands for com- 

li 

ponent i in the interval (0,T ) ; p^(x) = Prob[X^(T^.) = x] ; P (x) = 

Prob[X^(T .) _< x] ; and H^(x) = Prob [X^(T_^) _> x] . As indicated in the previous 
chapter, we assume that 

(A T ) X exp(-A T ) 

-L J- J. -L 

s) *7 

where is the total expected demand over ( 0 , T ^ ] , and is estimated 

from the technical replacement factor (TRF) and the installation schedule. 

MODEL I— UNITS SHORT 

The first model considers the objective function to be minimization of the 
expected number of stockouts over the interval ( 0 , T ] . Thus, we want to find 
those integer values s.,s 9 ,...,s which minimize the total expected units 
short subject to a constraint on the total procurement funds. If s. units 
of component i are stocked and x. units are demanded, then the total 
number of units short for component i is 



13 



x. - s. if x. > s. 
11 11 



0 otherwise. 



The expected number of units short for component i is therefore 



[ (x i ■ s P p-; (x ) » 

x . =s .+1 
i i 



and the optimization problem is: 



n °° 

Minimize Z(S) = 7 7 (x.-s.) p.(x.), (5.1) 

. i L i 1 1 1 

i=l x . = s .+1 

l l 

n 

subject to £ c.s. £ B, 

i=l 11 



and 



s . 

l 



nonnegative integer. 



This constrained nonlinear separable optimization problem can be solved 
using the generalized Lagrange multiplier approach. Let S be the n-vector 
(s^,S 2 » . . . , s^) , and let L(5;9) be the Lagrangian function; 



L(S;9) = 



n 

7 



n 



V 



/ / (x. - s.)p. (x.) - y(B - ) c.s ) 

■ 1 1111 . T 1 1 

1=1 X . = S .+1 1=1 

1 1 



Because of the separability, this can be rewritten as 



n 

L (S ; 0) = l L.(sJ - OB; 

i=l 1 1 



where 



19 



L . (s . ) = I (x.-s.)p.(x.)+0c.s., 

11 L , , 1 i r ii 11 

X .=s .+1 
1 1 



If the demand distributions were continuous, the optimal solution to (5.1) 
would be given by solving the n +1 simultaneous equations for s*, i = 1 
to n, and 0 * which satisfy [ 12 ]: 



3L . 

l 

9s . 
i 



= C ; i = 1,2, . . . ,n, 



9L 

90 



0. 



Because of the integer nature of the decision variables, the above procedure 
has to be modified slightly. Let 

AL. (s.) = L. (s.) - L. (s . - 1) 

ii ii ii 

be the change in the function L. resulting from an increase in the stockage 
level for component i from s_.-l t0 The optimal value of s_^ is 

the largest value of s^ such that AL^(s) < 0. Simplifying, we determine 
that 

AL . (s . ) = 0 c. - H.(s.) 

ii ill 

Therefore, s? is the largest nonnegative integer such that 

H.(s‘ f ) > 6 * 0 . i = 1,2, ...,n. (5.2) 

ii i 

The above simultaneous equations cannot be solved in closed form for the 
decision variables. However, tne solution can be obtained easily using the 
approach by Everett (Reference f 8 ]). Select a trial value for 0 and find 
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s*(9) from equation (5.2) for i = l,2,...,n. Then evaluate 

n 

B(0) = l c.s*(0) 

i-1 1 1 

If B(0) is close to B (within prespecified tolerance limits) the process 
stops and s_^ = s_^(0) . If B(0) is not close to B, a new value of 0 is 
selected and the process repeats. 

Aid in the selection of the new values of the multiplier is available 
from the theory of the generalized Lagrange multiplier approach presented by 
Everett [8]. The multiplier 0 has the mathematical interpretation of being 

a 3Z 

9 SB • 

Thus 0 can be thought of as a measure of the decrease in the total units short 
per additional budget dollar. Since the objective function is a monotone non- 
increasing function of B, 0 must always be nonnegative. If the amount of 

a current solution does not use up all of B then decrease 0 to increase 

the budget expenditures; if it exceeds the budget then increase 0 to de- 
crease the expenditures. 

The theory also provides guidance at each iteration as to the quality of 

the incumbent solution. A lower bound on the optimal value of the objective 

function is provided (Daeschner [5 ]) at each iteration from the inequality 

Z(S*) > Z(S*(0.)) - 0.(B-B(9.)) 

1 J 3 

Comparison of the present value of Z; namely, Z(S*($)), with the greatest 
lower bound yet experienced (for cases where the budget constraint was satis- 
fied; i.e., the amount used is B(6.) £ 3) will provide information as to 

j 

the maximum benefit to be gained by continued searching. For example, suppose 
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that the first three iterations yielded values of the above lower bound of 
(1003, 1155, and 1167). Also, suppose that the current value of the objec- 
tive function is Z(S*(9^.)) = 1168. We know that the optimal solution can- 
not be smaller than 1167, so there is little to be gained by continued 
searching . 

The iterative process described above will determine the values of the 
decision variables which will minimize the total expected units short subject 
to the given budget constraint. Exact equality of budget expenditures may be 
impossible to obtain because of the integer nature of the variables. There- 
fore, the result determined may not actually be optimal, but the difference 
is not likely to be significant. It is important to appreciate that, for a 
given budget expenditure B(6), the solution obtained is optimal. 

MODEL II— MODIFIED UNITS SHORT 

The model of the previous section minimized the total number of expected 
units short for the provisioned system subject to the budget constraint. It is 
a trivial matter to modify that model to allow for weighting the components 
by essentialities. Thus, consider the objective function 



Z 0 (S) 



n 00 

I l 

i=l x .=s .+1 
i l 



E i ( x i -s i ) p.(x.) 



obtained by weighting the units short for component i by the essentiality 
code E_. where E ^ > E^ if component i is judged to be more essential than 
component j. The minimization of Z^(S) subject to the budget constraint 
gives a formula for finding the optimal solution which is similar to that 
found in the previous section. The optimal stockage level for item i is the 
largest nonnegative integer s^ such that 
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Similarly, the units short objective function can be modified to mini- 
mize the total expected requisitions short. Let g_^ be the expected quantity 
demanded per requisition for component i, and consider the objective 
function 

n co ( x - s . ) 

Z 3 (S) = I l - 1 - p (x )) 

i=l x .=s .+1 ^i 

1 1 

obtained by dividing the expected units short by the expected requisition 
size to yield the expected requisitions short. For this problem, s* is the 
largest nonnegative integer such that 

i l i°i 

The procedure used to search for the value of 0 for these two modifications 
is the same as that discussed for Model I. 

MODEL III— SMA 

To conclude this chapter on units-short models we next consider an objec- 
tive function related to those described above, but which is stated in terms 
more closely related to presently reported measures of supply effectiveness; 
namely, supply material availability (SMA) . 

SMA is defined as the fraction of all requisitions that are satisfied 
by stock on hand over a given time period. It would seem reasonable to at- 
tempt to maximize the predicted SMA for the items in a provisioning package 
subject to a constraint on the total procurement investment. 

Let D_. be the total demand for item i during the provisioning 
interval (0,Tj . We assume that requisitions are for :>ne unit. Let PL be 



the number of demands satisfied, and let IK be the number of demands 
unfilled. Then, clearly, 

D. = F. + U. 
a 11 

and the expected number of demands satisfied (fills) is 
E (F ± ) = E(D ± ) - E(U ± ) = A i T i - E(U i ) 

Summing over the n items in the provisioning package, we get the total 
expected number of ’'fills’ 1 to be 



n n n 

y e(f . > = y a.t.- y ecu.) 

• S 1 i i i 

i=l i=l i=l 



Division by 
(SMA) . 



n 

A.T. yields the ratio of expected fills to expected demands 

i=l 



SMA 



n 

y E(F.) 
i=l 1 





00 

l 

X . =S .+1 
1 1 



I 

i=l 



(x . - s .) p. (x.) 
1 111 



A.T. 

i i 



( 5 . 3 ) 



Observe in (5.3) that the denominator does not depend on the decision variables 

and the summation term in the numerator in the total expected number of units 

short. Thus, the allocation of a given budget which minimizes expected units 

short will be the same as the allocation which maximizes SMA. However, since 

the value of SMA is likely to be more operationally meaningful than the total 
expected units short, the performance prediction should probably be stated 
in terms of SMA, This will aid in the communication of model results. 
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Just as with the units short model, the SMA model can be modified to 
include essentiality weights. All that is required is to weight the units 
short by the essentiality values and to modify the denominator to include 
also the essentiality values. The essentiality-weighted SMA (ESMA) would 
then be: 



ESMA 



i=l 



E. 

1 



x.=s .+1 
1 l 



(x.-s.^iCx.) 



n 

) 



i=l 



E.A .T. 

ill 



This overall value can also be written in terms of the ESMA's for the individual 
items. Let SMA^Cs^ be the expected SMA value for item it when s units 
are stocked. Then ESMA can be written as follows: 



ESMA 



n 

I 

i=l 



E.A.T . 

ill 



SMA_.Cs.) 



n 

v 



L 



i=l 



E.A.T. 

ill 



This alternative form illustrates that the overall ESMA value is a weighted 
average (weighted by essentialities and the expected number of demands) of 
the individual item SMA ' s . 

SUMMARY 

In this chapter we have considered models which allocate a given budget 
amount so as to optimize system performance with respect to units short or 
essentiality weighted units short. No consideration is given in these models 
to the length of time that must pass before a customer’s requisition can be 
satisfied when a stockout occurs. This time delay is uppermost in the minds 



Id 



of the customers as it can directly affect the readiness of weapon systems. 

The models developed in the next chapter remedy this shortcoming by explicitly 
considering time-weighed units short. 



6. TIME-WEIGHTED UNITS SHORT MODELS 



We develop two models in this chapter which consider not only the 
shortages which accumulate over (0,T], but also the length of time that 
each shortage exists. Classical inventory cost models recognize the effect 
that stockout delays have on system performance by imposing a stiff penalty 
cost for stockouts which grows as the stockout time increases. The models 
that we develop here do not assess penalty costs since they are virtually 
impossible to quantify. Instead, the objective function is the direct mini- 
mization of expected time-weighted units short (TWUS) . Such an objective 
function is appropriate when the readiness of the operational units suffers 
as long as replacement parts are unavailable. In a readiness measure such 
as operational availability, the time until a failed component is restored 
to serviceable condition is a key ingredient. When spares are available, that 
time is short. But when spares are not available, that time will be very 
long. 

These models assume that if there is a spare item in the system then 
it is '’immediately available” and that if there is no spare the demand be- 
comes a backorder which is filled by the next order which arrives from the 
manufacturer . 

MODEL IV— TWUS 

Let us restrict attention initially to a single component. Let X(T) 
be the random number of demands for the component in ( 0 , T ] . Let us assume 
that X(T) = m, and let 

0 < T-. ^ T 0 < ... < T < T 
12 m — 

be the times at which the m demands occur. If s units of cne item are 



stocked, the total time-weighted units short, given X(T) = m and 
T^,T 0 , . . . ,T m> can be described by 



I 



j=l 



(T - T 



s+j 



)• 



Figure 6.1 provides a graphic illustration of this derivation. The backorders 
are those demands giving negative net inventory. The length that each is 
outstanding is the time from when its demand occurs until the next order 
arrives at time T. 

NET 

INVENTORY 

s 



0 



m-s 

Figure 6.1: NET INVENTORY OVER TIME 

The conditional expected value of the time-weighted units short is therefore: 




E [TWUS (s) ' X (T) = m] 



0 



m < s 



m-s 

l (T-E[T s+j |X(T) =m]) if m > s 



( 6 . 1 ) 
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Let us now determine E[T^|x(T) = m] . Assuming a homogeneous Poisson process, 
the probability that the kth demand will occur between t and t+dt can 
be determined as follows: 



f_ ( 1 1 m) dt = Prob (t < T < t+dt X(T) = m) 

T k R ' 

P(k-1 events in (Q,t], 1 event in (t,t+dt], m-k events in (t,Tj) 

P (X (T ) = m) 



(At) k 1 exp (-At), (A(T-t)) m k exp(-A(T-t)) 

(k-1) ! (A dC) (m-k)! 

(AT) m exp(-AT) 

m! 



m! t,k-l ,T-t.m-k dt 

(k- 1) ! (m-k) ! ( T' ) 1 T ’ T ' 

r> 

Now, E[T ,X(T) - m] - / t f T (tjm)dc 

k 

('£ mj £ k-1 T-t . m-k dt 

J 0 C (k-1) ! ( m-k) ! V ^ T ' T 



The change of variable, u = t/T, simplifies the above expression tc 



T In “ 



m ! 



0 (k-1)! (m-k)! 



iC— 1 , \ m iC , 

u ( 1— u ) du 



The integral is recognized to be 
random variable with parameters 
conditional expectation that we 



T (— — ) 



the expected value of 
n = k and 6 = m-k+1 
seek is 

Tk 

m+1 



a beta distributed 
Consequently, the 



Returning to expression (6.1), we have determined the conditional expectation 



) ' 

I 0 




if m > s 



E[TWUS| X(T) = m] 



if m < s 



Expanding the summation gives (for the case m > s) 

T (m-s) (m-s+1) 

2(m+l) 

Finally, on summing over all possible values of m, the total expected time- 
weighted units short for the component is found to be 



Now that we have derived an expression for the expected time-weighted 
units short as a function of the stockage level s, let us write the optimi- 



weighted units short for item i when s. units are provisioned. We 
want to minimize the total essentiality weighted, time-weighted units short 
subject to a budget constraint. The Lagrangian statement of the problem is: 



oo 



E [TWUS (s) ] = l E[TWUS(s) |N(T) = m] Prob [N(T) = m] 

m=0 




(AT) H(s) - 2sH(s+l) + S - ( - S y 1 . - > H(s+2)} 



|{H(s+l) [At - 2s + S - ( - g 1 - ) ] p(s)(/\T-s)} . 



zation oroblem that we want to solve. Let V.(s.) be the expected time- 

1 i ^ 



n 

r 

) 



n 



Minimize L(S;0) 



i=l 



E.V. (s.) - e[3 - y c.s. ] . 

ill . l l 



i=l 
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^ et = ^i V i^ S i^ + ^ C i S i* before, the optimal solution 

s* is found by determining the largest s. such that AL.(s.;0) < 0 where 
n 1 11 

0 is such that £ c s. = B. After simplification AL.(s.;0) can be 
i=l 1 1 11 

written as: 



s . 

AL i (s.;0) = E.[(l-P(s i ))(T i - ^) + T.p(s i )] - Sc. 



Thus, s~ is the largest s. such that 



(l-PjU ))(T i - ii) + T p (s ) < ti. 

i i 



We solve for the optimal values using exactly the same search (on 9) 
procedure as discussed in the previous chapter. 

MO DEL V — MSRT 

One of the most widely reported measures of supply performance is the 
mean supply response time (MSRT) . This measure inherits its popularity at 
least partly through the role it plays in the expression commonly used to 
define availability (see Chapter 7) . In addition, MSRT is important by it- 
self as an indicator of the success of the supply system in meeting response 
time goals. As a measure, it considers both the likelihood of satisfying 
demands from stock on hand and the length of the delay in satisfying demands 
when the system runs out of stock. 

Let us determine the allocation of a fixed budget which will minimize 

the MSRT averaged over all of the items in a provisioning package. Let t 

be the time at which a failure occurs in (C,L], and let Y^(t) be the 

total number of failures which have occurred from FOC up to, but aoc ^ncluain 

t. Let R.(t) be the supply response time for the aemand occurring at time 
i 

t. We assume that R.^t) is a constant k. (possiblv zero) if Y.(t) s . - 

i l i i 



If the on-hand stock is not positive when the demand occurs at time t, then 
satisfaction of the demand will be delayed until units of time beyond 

receipt of the first replenishment order which we assume will occur at time 
T, Thus, 



( "i 

;<(t) - 



if Y (t) < s ± 



k. + (T-t) if Y (t) > s . 
11 i — 1 



The conditional mean supply response given the failure occurs at time t is 



S[R.(t)] 



k i ? [ Y i ( t ) < s i ] + (k + (T-t))P[Y i (t) >_s ± ] 

k. + (T- t)P [Y . ( t) > s . ] 

1 i 1 1 



The unconditional mean supply response time is found by integrating the above 
conditional expectation over (0,1] with respect to the density of the time 
t , But, it is well known (see Ross [21]) that the time of failure of an 
arbitrary event generated by a Poisson process is uniformly distributed on 
the interval ( 0 , T ] . Thus, the density of t is 



r\(x) 

L 




0 < x < T. 

— — i 

otherwise 



and the mean supply response time is 



MSRT.fs.) 
i i 



i . 

= / (k i + (X-tpH.Cs.)) I dt 



T. 

l 



= k. + I / (X-t.)H. (s.)dt 

i T. J n ii l i 



i 0 

I . s . (s .+1) s . 

= k. . (s . ' - — H. (s.+l) + — H . ( s . +2 ) 

•< 2 ' i \ t i ^2 ll 

“ " i " 2A . T . 

i i 
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Recalling the expression for the expected time— weighted units short we 
observe that 



MSRT^s ) 



k. + 

i 



V S 1 > 

A .T . 

1 i 



where V^(s_^) is the expected time-weighted units short when s_^ units of 
item i are stocked. Finally, the aggregate mean supply response time for 
all items in the provisioning package is obtained by weighting the individual 
means by the expected cumulative demand: * 



n 

y ( A . T) MS RT . ( s . ) 

i=l 1 * 11 

MSRT = — 

n 

" A.T . 

. - , li 



n n 

y k.A.T.-r y V . (s . ) 
. u , 111 . L n 1 1 
1=1 1=1 



l A.T . 
i-1 1 1 



(b . 2) 



Since the denominator and the first summation in the numerator do not 
depend on the decision variables, s^ through s^, it is clear that the 
stockage levels which minimize the aggregate MSRT are the same as those that 
minimize the total expected time-weighted units short. As with the TWUS model, 
the MSRT model can be modified easily to incorporate essentiality weights. 

All that is needed is to replace the weights 3 _^T_. in (6.2) by E^A^.T^ . 

Even though the TOTS model and the MSRT model will yield the same solu- 
tions, MSRT appears to be more operationally meaningful than is time-weighted 
units short. The next chapter develops a model which incorporates MSRT to 
yield another measure of performance of interest to the operational community; 
namely, operational availability. 
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7. AVAILABILITY MODEL 



The final model which we develop as a candidate for provisioning is an 
availability model which incorporates as a key ingredient the MSRT discussed 
in the previous section. We believe that the Navy Supply System would be 
better served by using one of the previous models, but we include the availa- 
bility model here to satisfy the present clamor for models which "spare to 
availability" . 

The commonly used definition of the availability for a single component 
is 



A . = MTBF . / (MTBF . + MTTR . 4- MSRT . ) 

l 1111 

where MTBF. is the mean time between failures, MTTR. is the mean time re- 
l i 

quired to repair /replace a failed unit when a spare is available, and MSRT_^ 
is the mean supply response time (see References [3 ], [14], [15], [23]). 
Although this expression is actually a limiting result which is valid only 
under certain conditions ([13], [21]), it is useful for providing insights 
into the factors which affect availability and of the sensitivity of availa- 
bility to those factors. For purposes of provisioning, we write the mean 
supply response time as a function of the total number of units provisioned, 
MSRT^(s^), as was done in the previous section. The availability expression 
then becomes 

A . (s . ) = MTBF./ (MTBF. + MTTR. + MSRT.(s.)) (7.1) 

li ii i ii 

Clearly, a primary objective of the Supply System should be to allocate 
resources so that weapon system availability is maximized. The key phrase 
^n the previous sentence is "weapon system." This is where one of the major 



difficulties arises in trying to develop availability models The availability 
of a weapon system is a very complicated function which depends on 1) the 
availabilities of its components, 2) the configuration of the weapon system 
(reliability block diagram), and 3) the way that the system is deployed 
operationally. Without making some simplifying assumptions, it would probably 
not be possible to develop any useful general purpose algorithms for maxi- 
mizing weapon system availability. 

Most existing availability models (see References [3] and [23], for 
example) make the following assumptions: 

1) Component failures are generated by a Poisson process (lifetimes 
are exponential). This means that MTBF^. = 1/A^. 

2) Component availabilities are given by expression (7.1). 

3) The weapon system is configured with the components all in series. 

4) The failure of one component is independent of the failures of 
the others. 

These assumptions lead to the formula below for the availability of the 
weapon system. 



a 

A (s ,s ) = n A. (s . ) (7.2) 

old n . , l l 

i=l 

The accuracy of the above expression depends on Che extent to which the 
assumptions are valid. We doubt that they are all satisfied for any real 
weapon system, and we doubt that the resulting estimated weapon system availa- 
bility would accurately reflect the actual observed system availability. As 
a consequence, we will refer to this model as the "pseudo-availability model." 
Nevertheless, as a model of system availability, ( 7 . 2 ) may be useful for 
the purpose of deciding the range and depth of items to be provisioned. 
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Consider the problem: 



maximize A 
s l’ * • • ,s n 



n 

n MTBF . / (MT3F . + MTTR. +MSRT.(s.)) 
. , 1 i 1 li 

1=1 



(7.3) 



n 

s. t . )* c . s . < B 

i=i 1 1 ~ 



= 0 , 1 , 2 , . . . 



Taking the logarithm of both sides of the objective function, we get: 



n 



in 



A(s- , , . . • , s ) = I £n(MTBF./ (MTBF. + MTTR. + MSRT.(sJ)). 

1 z’ n ii i i i 



i=l 



Let L(S.;9) = in (MTBF . /MTBF . + MTTR .+ MS RT .( s .) ) - 9c. s.. Solving problem 

l ii l ii ii 

(7.3) is equivalent to solving the Lagrange problem: 



max L(s,,...,s ;0) = y L.(s.;9) + 6B. 

1 ’ n’ i i’ 

i=l 



The optimal solutions are given by determining the largest values of s_^ 



suca tnat 



AL . ( s . ; 9 ) = L . (s . ; 9) - L.(s.-1;8) > 0. 

ii ii ii 



Simplification shows that s~ is the largest non-negative integer such that 



J. + MSRT.(s.-l) 

_ j_ i l 

J . + MS RT . ( s . ) 

l ii 



exp(9*c^) , 



where J. = MTBF.+MTTR. and 9* is that nonnegative Lagrange multiplier 
n 

value such that y c.s* = B. 

i=i 1 1 

Unlike the previous models, it would not make sense to weight the indi- 
vidual components by their essentialities in the availability model. In 
fact, the. series-system assumption is tantamount to an assumption of equal 



36 



essentiality of components; i.e., the failure of any component will render 
the system unavailable. True availability models which explicitly consider 
system configuration automatically handle essentiality by consideration of 
the impact of each component failure on system performance (see Jee [13]). 

No external weighting or assignment of essentiality codes would be necessary. 
Such "true" availability models are, however, very difficult to develop and 
implement . 

Because of the assumptions made in developing the pseudo-availability 
model in this section, it is unlikely that the system availability predic- 
tions provided by the model will be correct. Indeed, for even relatively 
small systems, the predicted availability will necessarily be very small. 

For example, a system composed of 25 components, each with an availability 

? 5 

of 0.95, will have a predicted availability of only (.95)" = 0.2774. 

Nevertheless, the budget allocation decisions made by the model may be very 
good. 

To conclude this discussion about the pseudo-availability model we would 
like to comment on a claim frequently made in the availability literature 
(see, for example, [ 3] and [23]). Some of the availability analysts argue 
that when sparing to availability, the only term in the availability expres- 
sion which involves the number of spares is MSRT . Therefore, they claim, 
maximization of system availability with respect to the number of spares 
for the components is equivalent to minimization of MSRT. A simple counter- 
example given in Appendix C shows that the claim is not true! However, for 
realistic choices of the various parameters, the allocations determined 

from the two models ("MSRT and A ) may be very much alike. The numerical 

o 

evaluations which follow in the second report examine this issue. 
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8. COMMENTS ON THE MODELS 



INTRODUCTION 

The previous chapters have presented three basically different provision- 
ing models: 

(1) units short (SMA) , 

(2) time-weighted units short (MSRT) , and 

(3) pseudo- avail ability. 

In this chapter we provide brief comments about the models as related 
to three areas: 

(1) incorporation of a random protection interval, 

(2) theoretical properties of the models, and 

(3) solution of a related problem for use as a budgeting tool. 

RANDOM PROTECTION INTERVAL 

The previous models ail assume that the length of the protection inter- 
val (?0C to POC-HT] is known. In reality, T is a random variable since both 
the procurement leadtime and the time from POC to the first replenishment buy 
are random variables. If one had information about the distribution of the 
random variable T, the previous models could be extended (at the cost of 
additional mathematical complexity) to accommodate the randomness of T. 

Let ( t ) be the probability density function for T. What is required 
in the previous models to incorporate the randomness of T is to replace T 
everywhere by a dummy variable, say u, and to integrate, weighting by the 
density g T (u), over all possible values of u. One common assumption for 
the form of g^(u) is the gamma density. The gamma probability family has 

j. 

the proper range of values (0 to °°) , and the family is extremely rich 
in the shapes that it can assume. Many theoretical and empirical studies 
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have shown that the gamma distribution frequently approximates well the 
empirical distributions of leadtime. Furthermore, when the gamma distribu- 
tion for leadtime is combined with the Poisson distribution for demand in an 
interval, the resulting distribution of leadtime demand is negative binomial 
(see Chapter 3 of [12]). This theoretical fact accounts for the selection 
of the negative binomial as one of three distributions used by the Navy ICP T s 
in some of its wholesale provisioning and replenishment models. 

For purposes of implementation in the provisioning models it is unlikely 
that any information will be available to adequately estimate the distribution 
of T. Also, since inclusion of a probability distribution for T increases 
the complexity of the models, we suggest that T continue to be treated as 
a constant. Of course, every effort should be made to estimate T as well 
as possible. 

One of the elements of T is the length of time from POC to the time 
that the first replenishment buy is made. In addition to the elements already 
included in the models, this time period depends on an external element not 
yet considered — the replenishment reorder level. Given the size of the first 
provisioning buy, say s, and the size of the first replenishment reorder 
level r, it is easy to determine that the time from POC to the first replenis 
ment buy is Erlang distributed (assuming a homogeneous Poisson process with 
rate A generates demands) with parameters n = s-r and A. This is because 
the delay is the amount of time required to experience s-r demands with 
exponentially distributed times between demands. The problem is that in order 
to integrate the provisioning and replenishment models it would be necessary 
to include the replenishment reorder level r as yet another variable in the 
provisioning model so that the size of the provisioning buy depends on r. 
Furthermore, the value of r in such an integrated model would logically 
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depend on the provisioning stockage level so that solution of the model would 
require an iterative technique greatly complicating the model and increasing 
the computational burden. Thus, we suggest that the first replenishment 
buy not be triggered by the inventory position dropping below a reorder 
point, but be made automatically at a fixed time after POC, say POC+At. 

If the estimates of demand, installation schedules, and procurement leadtimes 
are accurate, it is reasonable that the first replenishment buy be made at POC 
(i.e.. At =0). Delaying the replenishment buy to some later date to allow 
the initial TRF estimates to be updated by actual demand data only lengthens 
the provisioning protection interval and increases the size of the required 
provisioning buy. 

THEORETICAL PROPERTIES 

The models developed in this report can be expected to lead to different 
solutions to the optimal budget allocation problem. In contrast to the Varia- 
ble Threshold Model, each of the models also provides a prediction of the 
expected performance to be achieved from a given solution. Furthermore, some 
theoretical properties of the models can be examined to determine which model 
best reflects our intuitive notion about how the system is degraded by the 
presence of shortages, about the effect of uncertainty in our estimates of 
the demand rates, and about the propensity to stock items as a function of 
their demand rates, unit costs, essentialities, and leadtimes. 

All of the models are more likely to stock the low cost, high demand-rate 
items, all other things the same. Weighting by essentialities is the only 
way offered here to oppose that tendency. Also, the models are all piece- 
parts oriented. No explicit consideration is given in any of the models for 
r, weapon system effectiveness". However, essentiality weights can be used to 
reflect the importance of the individual items with respect to weapon system 



effectiveness. Proper assignment of essentiality weights should reflect 
what happens to the readiness of a weapon system if a shortage occurs in an 
individual item. 



A COMPANION PROBLEM AS A BUDGETING TOOL 



In the models which we have developed, we have assumed that a budget B 
was given at the outset. DODI4140.42 requires that the provisioning budget 
be determined by using the COSDIF model. As discussed in Appendix B, the 
COSDIF model first determines if an item is to be stocked by comparing the 
expected cost of stocking the item with the expected cost of not stocking the 
item. If the former cost is smaller, the item satisfies the range rule. 

The depth of a selected item is the expected demand during a procurement 
leadtime plus one quarter. The total "priced-out" value for all items satis- 
fying the range rule is the provisioning budget (excluding the budget allowed 
for insurance and NSO items) . 

This procedure for determining the budget for a provisioning budget is 
cost-based requiring cost estimates that may not reflect accurately real 
costs. Furthermore, the COSDIF procedure does not provide any estimate of 
the performance that might be expected from the resulting budget. The 
alternative models presented in this report can each be rewritten in a 
mathematically-related form which could be solved to provide a tool useful 
for determining a provisioning budget. 

We saw that each of the models described earlier could be written in 
the primal form: 



min(max) 



y f.(s.) 

i=i 1 ^ 



subject to 




( 8 . 1 ) 



The companion problem can be defined as follows: 



minimize ) c.s. 

11 

n 

subject to £ f.(s.) (>) G*. 

i=l 

where G* is a specified performance goal. (The direction of the inequality 
depends on whether (8.1) is a minimization or a maximization problem) . This 
companion problem can either be solved directly using a generalized Lagrange 
approach like that used to solve the primal problems, or the primal problem 
can be solved for a range of budget values and the smallest budget satisfying 
the performance goal be selected. 

Solution of this companion problem will provide planners with a budgeting 
tool through which they can justify a budget on the basis of the performance 
that can be attained. Thus it provides a mechanism for relating resources 
to readiness. 



9 . SUMMARY 



We have developed alternative wholesale provisioning models which are 
intended for use by both SPCC and ASO. The models determine how to allocate 
optimally a given provisioning budget with respect to various measures of 
performance. The measures of effectiveness can be summarized into three types 

1) Essentiality-weighted units/requisitions short or SMA, 

2) Essentiality-weighted time-weighted units short or MSRT , and 

3) Availability. 

In addition to providing the optimal budget allocations, each model provides 
a prediction of the performance to be expected during the provisioning period. 

The models all assume that a budget, as determined by the COSDIF process, 
is given. However, each model can be converted easily into a related optimiza 
tion problem which determines the minimum budget required to achieve a speci- 
fied performance goal. We feel that it makes more sense to determine a budget 
in this manner, than by using the COSDIF approach which is heavily cost 
oriented and which provides no estimates of the performance to be expected 
with a given budget. 

The alternative models solve simultaneously the range and the depth 
problems. Thus, no additional model is required for the budget allocation 
process . 

In order to evaluate the proposed alternative models, the models should 
be exercised with sample data, and the results examined carefully to observe 
how each model allocates the provisioning budget. The results of each model 
should be tested with respect to several measures of performance (units short, 
SMA. time-weighted units short, MSRT, availability, etc.;. One should then 
select that model which provides results that best agree wi tn a decision 



4 3 



maker’s perception of how the system should perforin. In exercising the 
alternative models the present models should also be tested to serve as a 
baseline for comparison. This is the type of model evaluation that will be 
discussed in the second report. 

However, even before a real data performance evaluation is performed, 
the alternative models can be evaluated theoretically or intuitively to deter- 
mine if they attempt to achieve the goals of the supply system. We have com- 
mented on various properties of each model in the sections developing the 
models and in the previous chapter. Because of those theoretical properties 
of the models, and because of the role that time-weighted shortages (MSRT) 
plays in the expression for availability, we believe that the MSRT model makes 
the most sense for provisioning. It captures the essence of the objectives 
of the supply system while avoiding some of the problems and data requirements 
that would be associated with the pseudo-availability model. 



APPENDIX A: INITIAL DEMAND FORECASTING 



THE PRESCRIBED PROCEDURE 

The initial demand prediction process has been delineated in Reference 
[7] and begins with a forecast of "program data'* from POC out over the two- 
year demand development period. This data is a month-by-month schedule of 
the anticipated end item installations over that time period. The time- 
weighted average months program (TWAMP) is computed next for the appropriate 
program time base (PTB) . Reference [7] allows the PTB to have a maximum 
length of one year. Shorter PTBs are 3 months and 6 months and are appro- 
priate for computing needs for expensive spares. However, SPCC and ASO use 
only a 12 months PTB regardless of an item's value because of the workload 
created by having shorter PTBs (a review is needed at the end of each PTB 
interval for the shorter ones) . 

The concept behind the TWAMP is quite simple. Consider Figure A.l; the 
TWAMP is computed by determining the area under the curve of the installation 
schedule and then dividing it by the length of the PTB. This results in an 
average number of end items to be supported over the PTB. This average is 
denoted as the initial TWAMP or TWAMP by Reference [19]. 

Jt 



J — 1 ! 1 > 



Figure A.l: 



An Installation Schedule 



The following formulas are provided by Reference [7] for computing the 
initial TWAMP. They obtain the area by summing the vertical slices of 
the area spanning each month. Any installations occurring in a specific 
month are assumed to begin operation in the middle of the month. 

( I i f 

I — for m = 1 

m ) m-1 I 

I l I, + for m >_ 2 

V k=l 

where 



I, = the number of units of the end item to be installed 
in month k. 

PTB 



l 



TWAMP , 



m=l 



m 



PTB 



The initial annual demand rate for a given spare part can be computed 
by using the following formula [19]: 



D, = TWAMP * N * BRF 
I 



where 



N = the number of units of a given replaceable part in 
the end item. 

BRP = tne ’’best replacement factor” which is, in fact, the 
estimated failure rate of one unit over a year. 

The BRF is initially based on a technical replacement factor (TRF) which 
is the contractor’s estimate of the attrition rate. 

The COSDIF formula requires an initial estimate of the steady state 
annual demand rate D^ c . This is computed by first determining the number 
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of end items to be installed over the first two years (its value is denoted 



as TWAMP^ in Reference [19]). Then 

D ss = TWAMP SS *N*BRF. 

If a spare part is repairable then the formula considers only attrition 
demand. In that case. 



D ss = TWAMP *N*BRF[1 - CR_R*RSR] (A.l) 

where : 

CRR = the carcass return rate, and 
RSR = the repair survival rate. 



The initial estimate of demand for repairables will be introduced 
later. 

The development of the provisioning budget constraint of Reference [7] 
requires that the initial demand during the procurement leadtime (PCLT) 
plus one quarter be computed. This is then the depth value to be used to 
develop the budget constraint. Reference [19] also needs the demand during 
PCLT for the Variable Threshold model. For consumable items the formulas 
are, respectively, 



PCLT 



= 0 . 



PCLT 



D 



PCLT4-1 



v- 



. (PCLT 4-1) 



For a repairable item the computation used in Reference [19] includes the 
repair turnaround time (TAT). The formulas are based on the following logic. 
Without any carcass returns the attrition demand is the same as Equation 



(A.l) with TWAMPgg replaced by TWAMP,. If the repair process is "up and 
running" at POC minus TAT, then the first repaired item will be available 
at POC and will provide a replacement for some unit of the item when it 
fails. The net attrition loss is then 



TWAMP *N*BRF*[1 - CRR*RSR] 

If, as is expected, the repair process can T t begin until some items fail, 
then the net attrition loss just described must be increased by the number 
of units required to "charge" the repair pipeline assuming that TAT is 
shorter than the PCLT or PCLT plus one quarter. The respective formulas 
for demand of repairables are therefore: 



PCLT 



PCTT TAT 

D,*[— r^*(i - CRRrtRSR) + — -*CRR*RSR] 

A a *4 



^PCLT-rl 



D *[ (?Cj ' T+1 L (1 - CRRaRSR) + ^LcRR*R5R] 

A 4 4. 
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APPENDIX B: THE COSDIF FORMULA 



B. 1 THE PHILOSOPHY 

The COSDIF formula was provided by Alan Kaplan, of the Army T s Inventory 
Research Office, who was a member of the team that developed the DODINST 
4140.42 [7]. It is derived from the Army 1 s wholesale range model. It can be 
viewed as the results of a simple decision model as illustrated in Figure B.l. 





r 

States of Nature | 


Decision 


No demand during DDP 


1 

Demand during DDP 


1 ! 




Make a 

provisioning 
i buy 

1 


Costs of procurement 
‘plus two years f holding 
'costs 
| 


Two years of average 
annual variable costs 

i 

I 


Do not make a 
provisioning 
buy 


No costs 

| 

1 


i 

Costs of spot buys 
during first year and 1 

average annaul variable 
costs for second year 

. . _ i 



Figure B.l: A Provisioning Decision Matrix 

The expected costs of each decision can be evaluated when all the costs 
and the probability values for the states of nature are known. That decision 
which corresponds to the least expected cost is then the optimal one. 

The COSDIF formula given below determines the difference between the 
expected costs of making a buy and not making a buy. 

If the COSDIF value is negative, then the costs associated with not making 
a buy are greater than those for making the buy, and is therefore optimal to 
make the buy. If the COSDIF is positive, then making a buy is more costly. 

If the COSDIF is zero, either decision would be optimal but not making a buy 
is less work for the provisicner so no buy is made. 
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COSDIF 



COSDIF = 


(D |D_) a {CP + 2*H*C*(R+D__/FP)} 
o T SS 

+ (1-D q |D t ) * {CP*FP + H*C*D s /( 2*FP) + D T *Cl} 
- (i-d o |d t ) * {d t *(csp + plt*a/4) + d t *c*p) 



wh er e 



D ss 


= steady state annual demand (attrition) 


d t 


= total steady state annual demand (TWAMP^*ERF) 


D o i D T 


= probability of no demand in two years, given an annual 
steady state demand forecast (total) of D . This 
value will be obtained from the conditional probabili- 
ties stored in the SCA after the computed is 

rounded up at .5. 


CP 


= cost of procurement [the small purchase value is used 
if the value of the buy quantity is less than or equal 
to the small/large purchase breakpoint; if not, the 
large purchase value is used. The specific values to 
be compared to the breakpoint are the values of the 
buy quantity (R + D„/FP) and the economic order 
quantity (D 55 /FP) for lines 1 and 2, respectively, 
of the formula] . 


H 


= holding cost rate 


C 


= unit price 



R = reorder point quantity 



FP 


= frequency of procurement 


Cl 


= cost of issuing stock 


CSP 


= cost of spot procurement 


?LT 


= production lead time 


S 


= implied shortage cost [the larger of the shortage cost 
for COSDIF (specified in the SCA cog constants) and 
the cost to hold, K*C] 


P 


= spot buy premium rate 


In examining 


che COSDIF formula we see that the second and third lines 



correspond to the outcomes of the alternative decisions when there is demand 
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during the DDP (see the last column of Figure B.l). However, they show only 
the average annual variable costs for one year and the expected costs of 
making spot buys for one year. This is because the average annual variable 
costs for the second year cancel each other when the cost difference is 
taken according to [7]. 

B . 2 A CRITIQUE 

The COSDIF formula serves as the range model in the development of the 
provisioning budget constraint. As a consequence it is appropriate to examine 
its elements and to critique the assumptions. The purpose of this critique 
is to address certain details of the formula so that the reader will be able 
to identify where improvements in the current use of the model by the Navy 
can be made. It is not possible to change the formula without approval, and 
recent attempts to do so have not been successful. These will be discussed 
below. 

A rather puzzling assumption involves the use of the steady-state demand. 

It is especially so when we realize that the depth model does not use the 
steady state value. Perhaps it is true that, because cost differences are 
being used, use of instead of should make little relative difference. 

This should be evaluated. 

The COSDIF model considers only two states of nature with respect to the 
occurrence of demands over the two year demand development period — zero demands 
or positive demands. SPCC uses probability estimates of these states which 
are a function of the expected mean demands. The probabilities are given in 
Reference [10] . The probabilities used by SPCC are more appropriate than those 
in [7] since the latter were based on inventory records for ncn-Navy items. 

In fact, [7] required each service to come up with better values within two 
years. The probabilities used by SPCC can be further Improved by developing 
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values for each two-digit cog as was attempted in [10] . This is strongly 
recommended since cog classes do represent items of a similar nature. It is 
interesting to note that ASO has apparently continued to use the probability 
values given in [7] for all of its cog classes [6], 

It could be argued that greater resolution in the demand distribution 
should be considered in the COSDIF expression — perhaps estimates of individual 
positive demand quantities for, say, quarterly or yearly periods over the 
two-year DDP . Even if the probability was kept at a single value, the form 
of the average annual variable costs to stock an item over the two years sug- 
gests that several buys will be made and orders received. The buys might be 
made but, as was noted in the main part of this report, no deliveries are ex- 
pected until after the end of the DDP. In fact, the average variable cost 
expressions assume that steady state is reached instantaneously at POC (as 
the use of D^^ implied earlier) . A switch from a no-stock to a stocking 
mode between the first and second year is virtually impossible because of the 
leadtimes. The spot buy situation will continue through ail of the DDP under 
current replenishment leadtime values. There will be no holding costs since 
no inventory will be on hand. 

The term of the COSDIF equation corresponding to the zero-demand state 
of nature is the closest to being correct ever the two year DDP. The impact 
of stocking and holding for two years is significant and is not unexpected in 
light of the concern in [7] about excessive unused stocks of provisioned items. 

Considering the spot buy term in more detail, it is interesting to note 
that there appears to be no explicit term for an item’s essentiality. In an 
earlier document (DODINST 4140.39 [6]), the replenishment model included such 
a term. If one searches through (7] for a reason why essentiality was not 



included, the only consideration that is found is in Paragraph IV. E. : 



Non-demand-based items, which do not meet the insurance item criteria 
for wholesale level stockage may be stocked in the wholesale system 
only if there is an overriding requirement to do so based upon their 
essentiality in a selected weapon system. In this instance, an item 
will be stocked as a Numeric Stockage Objective (NSO) item. 

Reference [7] also notes that 

no deviation to the basic formula is authorized, the use of cost data 
more relevant to a particular inventory should provide sufficient 
flexibility to account for the variability required in support of 
different systems. Data developed for ICP wholesale level systems 
in accordance with DODINST 4140.39 should be useful in establishing the 
most appropriate holding costs, procurement costs, and implied shortage 
cost . 

Additionally, a constraint is imposed such that the cost of a requisition on 
backorder must be at least equal to the cost to hold an equivalent amount of 
material. It may also be desirable to use a larger value for the holding cost 
in the COSDIF than is used in Reference [6] to reflect a higher obsolescence 
probability for marginally stocked items [7]. 

A variation in the COSDIF which includes an essentiality term was presented 
in Reference [9] . However, the goal of that reference was to justify the 
variable threshold model (which does include an essentiality term in the E.ISK 
equation used to determine the depth) as an alternative model for implementation 
of the budget constraint. It is important to note that D55, which is SPCC’s 
implementation of the COSDIF, does not use this variation in creating the 
budget constraint . 

In July 1981, the Office of the Chief of Naval Operations [4] proposed 
changing Reference [7] to correspond to the form of the COSDIF in [9]. The 
response from the Deputy Assistant Secretary of Defense (Logistics and 
Material Management) was [11] : 

While we understand the thrust of your initiative to revise DoD 
Instruction 4140.42, we do not agree with the necessity of making 
the requested changes. This instruction currently recognizes essen- 
tiality in the initial spares computation and allows stockage of essen- 
tial items even if they do not qualify for stockage through the use 
of the cost differential tables. 
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When a shortage cost is considered in any inventory model supporting a 
weapon system, it is implicitly a measure of the essentiality of the item. 

The value that can be used for A g in the COSDIF equation is not constrained 
by References [7] and [17] to be that of the replenishment model. In fact, 
the values used in the replenishment models are tied to a replenishment budget 
constraint in any given fiscal year. As a consequence, the shortage cost 
values in the replenishment model can change each year. It should be obvious 
then that such values are inappropriate for provisioning. 

Control over A^ will have a significant effect on the development of 
the budget constraint. The larger the A^ value, the larger will be the 
range of demand-based items used in developing the budget constraint. The 
bottom line is to obtain a larger budget constraint. The determination of 
better A^ values should be viewed as an important first step of improving 
the provisioning process. The ASO provisioning model attempts to do this. 

It adjusts A g in an effort to gain an SMA of 85% for the depth allowed by 
[7] for determining the budget constraint. This is done separately for 1R 
cog and for the group of 2R and 6R cog items. The ASO model develops a 
cost constraint for five different values of A^ using application/operation 
D54 [17] . For each A , the depths of the items are used to estimate the 

3 

SMA which would result. A curve is then constructed of SMA versus investment. 
Additional A values as needed are used to generate additional points as the 
process zeros in on the 85% SMA point. The investment value for the 85% SMA 
point is then established as the budget constraint. The actual range and 
depth of provisioned items are computed using D52 [18] . D52 uses a load list 

approach developed prior to Reference [7]; the goal of that approach is to 
minimize expected units short. The A^ value generated by D54 is used only 
to generate the budget constraint. D52 uses a Lagrange multiplier approach to 
implement that budget. 
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It is debatable whether SMA is the appropriate effectiveness goal to use 
in light of the current trend toward A . Additionally, SMA is measured only 
for those items passing the COSDIF range criterion. However, the approach 
used in the ASO model is certainly worth considering with availability goals 
used in place of the SMA goal. 









APPENDIX C 



This appendix gives a simple counter-example to show that maximization 
of ”pseuao-avaiiability H as defined by the relationship 

n MTBF. 

A = n 

o . MTBF. + MTTR. + MSRT ,( 3 .) 

11 i 11 

does not yield the same solution to the stockage allocation problem as 
minimization of mean supply response time (MSRT) or time-weighted units short 
(TWUS) . Thus, the claim that ’’sparing to availability” is equivalent to 
sparing to minimize logistics delay time is not true. 

Example: Consider a system consisting of two components having the parameter 

values given in Table Cl: 



PARAMETER 


1 


2 


MTBF 


0.200 


0.100 


PCLT 


1.0 


1.0 


MTTR 


0.0822 


0.0274 


C 


5. 


10. 



TABLE Cl: SYSTEM PARAMETERS 

where PCLT is the procurement leadtime in years; C is the unit cost; 

MTBF = l/X is the mean time between failures; and MTTR is the mean time 
to repair. Suppose we are allowed a total provisioning budget of B = S20. 
Then it is clear chat there are only 3 undominated feasible solutions to the 
allocation proolem: (4,0), (2,1), and (0,2). Table C2 gives the mean 



supply response times for each of the feasible values of stockage levels for 
components one and two. 



MSRT(s) 



s 

Item 


0 


1 


2 


3 


4 


i 


182.50 


124.00 


79.51 


47.80 


26.83 


2 


182.50 


149 .65 


120.45 


X 


X 



TABLE C2: MS RT vs. s 



From the information in Tables Cl and C2 it is easy to determine the overall 
mean supply response time and the system availability for each of the three 
candidate solutions. These values are given in Table C3. 



CANDIDATE SOLUTION 


MSRT 


A 

o 


(4,0) 


130.61 


0.0896* 


(2,1) 


126.27* 


0.0744 


(0,2) 


141.13 


0.0559 



TABLE C3 : SYSTEM MSRT AND A 

o 



The values in Table C3 reveal that the solution (2,1) is optimal with 
respect to total MSRT whereas the solution (4,0) is optimal with respect 
to system availability! 
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